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Exact chiral invariance at finite density on the lattice 
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We propose a lattice action for the overlap Dirac matrix with nonzero chemical potential which 
is shown to preserve the chiral invariance on the lattice exactly. We further demonstrate it to 
arise from the Domain wall by letting the chemical potential count only the physically relevant wall 
modes. 
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I. INTRODUCTION 



Our world of strongly interacting particles has two light quarks with masses much smaller than Aqcd, the scale of 
Q ' the corresponding theory. Quantum Chromo Dynamics (QCD). While the mass of the strange quark is comparable 
to this scale, other quarks are heavy. It has been argued on the basis of the corresponding symmetries of the theory, 
called the chiral symmetries, that QCD may have a critical point in the temperature (T)-baryon number density (or 
' equivalently, the baryonic chemical potential /is). Theoretical [H-Q as well as experimental searches for locating it 
(N , are currently going on [4j. Its discovery would be exciting in many ways. Apart from becoming a new milestone in 
our understanding of the nature of the strongly interacting matter, it would also be unique compared to the other 
known phase diagrams in the way theory and experiment compliment in locating the critical point in it. 

Due to the essentially non-perturbative nature of the problem, theoretical efforts based on a first principles approach 
employ lattice QCD techniques which have made successful predictions for many hadronic quantities such as the heavy 
jjj ■• meson decay constants. Indeed, it has also been very successful in application at finite temperature, leading to several 
f~| ' interesting results for the RHIC and other heavy ion experiments [Sj] • But its foray in the finite density domain has 
been hampered by serious conceptual problems. In spite of the so-called 'fermion sign problem', which refers to the 
complex measure of the theory for nonzero chemical potential, many attempts have been made to explore nonzero 
fiB, some claiming success [Ull, @| in locating it, while others arguing for a lack Q of a critical point. Most of these 
' computations use the staggered quarks, primarily since the critical point is thought to be related to chiral symmetry, 
■rj" Staggered quarks have a U{l)v x U{1)a chiral symmetry on the lattice, and a corresponding order parameter, the 
0^ chiral condensate. Studying its variation with T or /is, one can look for a chiral phase transition. The critical point 

, would then just be an end point of a line of first order phase transitions. 
T-H ' Staggered quarks, however, i) break the flavour and spin symmetry on the lattice, ii) have no flavour singlet axial 
T-H I symmetry and iii) have a stronger [8] 'rooting problem' at nonzero /is. Since the QCD critical point is expected [9| 
. to exist if only if one has two light flavours, and the flavour singlet anomaly is mildly temperature dependent [lO|, 
" ^ I ' it appears desirable to improve upon them. The overlap Dirac fermions or the closely related domain wall 
J> , fermions [l^, offer such a possibility to improve. Indeed, the overlap quarks have all the symmetries of the continuum 
QCD and also have a an index theorem [l^ as well, raising the hope that even the anomaly effects could be well 
^ , treated. Unfortunately, adding the chemical potential turns out to be nontrivial for them. Bloch and Wett ig I14| 
' made a proposal to do so but it violates [lH the exact chiral invariance on lattice as does the simple addition jl6j| of 
a baryon number term. In this letter we propose an alternative which does have exact chiral invariance on lattice for 
any value of the lattice spacing and any chemical potential. It therefore can permit, in principle, the task of mapping 
out the T-^B phase diagram, assuming that the algorithmic developments can handle the fermion sign problem well. 
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II. FORMALISM 

The massless continuum QCD action can be written in a form where the chiral symmetry is manifest in terms of 
the fields appearing in the action: 



Sqcd = J d^x dri^pij - F'^-'F^JA] 



d\ dT[ V {i>^p^,) - F^^F^./A] , (1) 



i=LM 



where V'l = (1 — 75)'/'/2 and -0/?. = (1 + 75)V-'/2 with — ^'(l + 75)/2 and ipR = ■0(1 — 75)/2. The second term is 
the action for gluons which will play no role in our discussion below. We assume below that some usual convenient 
form has been chosen. Adding the canonical to the first line for investigating finite density effects is the same as 
adding ^•tpij'^ipi to the second line, leaving the manifest chiral symmetry intact. We propose that addition of chemical 
potential on the lattice be done in this explicit chiral symmetry preserving manner as well, opting therefore for the 
overlap quarks. 

The overlap quarks have all the symmetries of the continuum QCD but also have a nonlocal action: 

S'_F = ^ V'n aF>ov,nm 0m , (2) 

where the sum over n and m runs over all the space-time lattice sites, a is the lattice spacing, and the overlap Dirac 
matrix Dqv is defined by aDov — 1 + 75S5"-(75^w). sgn denotes the sign function. Dw is the standard Wilson-Dirac 
matrix on the lattice but with a negative mass term Af € (0, 2): 



1=1 



The overlap Dirac matrix satisfies Ginsparg- Wilson relation [17| . {75. D} = aD^^D and has exact chiral symmetry 
on lattice. The corresponding infinitesimal chiral transformations [l8| are 

Sij) = ia^l5{l-'^Dov)'4^ and (Ji/i = ia'0(l - |l'oi,)75 . (4) 

An alternate set of transformations, differing by terms of 0{a) is 

(50 = 1075(1 — aDoy)ip and StJj = iaipj^ . (5) 

Since one needs to write chiral projectors on the lattice to mimic eq.(IT|), we focus on eq.(IS]). The generators of the 
transformation in eq.® satisfy 7I = 1 and 7I = [75(1 — aDov)]^ = 1- On the other hand, the corresponding ones in 
eq.(|4]) do not square to unity. One therefore defines the left-right projections for quark fields as -0^ = (1 — 75)0/2 
and ijjfi = (1 -|- 75)^/^/2, leaving the antiquark field decomposition as in the continuum above. Such a decomposition 
is commonly done for writing chiral gauge theories (19i] on the lattice and is possible since and ip are independent 
fields in the Euclidean field theory. In analogy with the eq.([T|) of continuum QCD, the action for the overlap quarks 
in presence of nonzero chemical potential may now be written down as 



ov 

+ afij )-0„,i] (6) 

n i=L,R 

= Vi„[(l - a^i'y'^/2)aDoy + a^7'^]0„ . (7) 

n 

It is easy to verify that i) this action is invariant under the chiral transformation eq.(I5]) for all values of a/i and a 
and ii) it reproduces the continuum action in the limit of a with /i l^/M scaling. In order to obtain the order 
parameter for checking if the symmetry is spontaneously broken, one usually adds a linear breaking term. Adding a 
quark mass term, am^tpfitpL + '4'l'4'F!.): on.e obtains the order parameter valid for all T and fj, on the lattice by taking 
a derivative of the log of the partition function with respect to am as, 

/0-0\ = lim lim (Tr - — — ; — ^ -7/7-^ — r^) ■ (8) 
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The only real eigenvalues of aDov are or 2, with only the former contributing to the order parameter. Defining 
a matrix Kg^ = -Dod(1 ~ o,Doy/2)~^ , such that {'y5,Koy} = 0, the order parameter can be written in a form more 
analogous to the continuum : 

i^^) = (Tr ^ i) . (9) 

Although the discussion above is for a single flavour of quark, i.e, U{1)l x symmetry, its generalization to 

Nf flavours is straightforward. Indeed, since it relies only on the spin-structure, the flavour index as well as the 
corresponding generator matrices just carry through. 

III. DOMAIN WALL FERMIONS 

The action in eq.Q was obtained by demanding the left- right symmetry to be explicit. A physically more intuitive 
way to arrive it is in the domain wall formalism. Since only the massless domain wall modes are physical although 
there are many massive unphysical modes, the appropriate way to introduce is as a Lagrange multiplier for the 
number of these massless modes. The domain wall action of [20| then for nonzero chemical potential, /x, is 

^-^ ^-^ L a 

x^x' s,s' — 1 

-t- an 74 5^^^, {Ss,iSs',iP- + P+Ss,NrJs'.NrJ + am Sx,x' (S^^iSs' ,N5P+ + PSs^NrJs' .i)] fpix' , s') , 

where P± — (1 ±75)72 and and as are the number of sites and the lattice spacing in the fifth direction respectively. 
The physically relevant 4D massless fermion field is identified with the fermion fields at the boundaries of the fifth 
dimension as, 

i^^P.lPl+P+i^N, , ^^i^lP++i^N,P- ■ (11) 

It is convenient to visualize the five dimensional action in terms of the fields 7]i localized on four dimensional branes 
existing at each site i along the fifth dimension, as in [2l|. Defining a transfer matrix, T = {1 + a^Hw P+)~^ — 
a^HwP-), between pairs of neighboring branes, where Hw = ^bDw, the five dimensional action can be rewritten in 
terms of these fields as, 



S = 



1^771 - maP^ + afi {a^HwP^ - 1) ^ 74^^-^ Vi ^ Vn^T ^ — maP^ ~ afj, {a5HwP+ + 1) ^ 74^+) Vi 



N5 



(12) 



The fermion fields rji can be integrated out successively, resulting in a partition function of the form, 

Z5D = J VUe-^'=' J detD(5)(ma,aAt) , (13) 

where J is the Jacobian for the transformation from tp-fi.elds to yy-fields and the five dimensional determinant D^^^ is 
given by 

2?(^) = det F_ - maP+ + afi {a^HwP- - 1)"^ JaP- - (^P+ - maP^ - a^i {ar^HwP+ + 1)"^) liP+ ■ (14) 

In order to obtain the overlap matrix, the contribution of the bulk five dimensional modes needs to be removed from 
the partition function. Following [22[, we introduce pseudo-fermions and obtain the partition function of interest as 

„ X /"^T. det D^^^ma, au) , , 

ZQCD{T,fi) = / VU e-s^ j-f^T^w ^ ; ' 15 

J det D(^^^(ma = 1, a^ = 0) 

where detD^^\ ma = 1, a/i = 0) is the contribution from the pseudo-fermions. Taking first as — > limit and then 
N5 — >■ 00, the ratios of determinants turns out to be det[Do« + (1 — Dov/'i){ma -\- a/x7^)], where the dimensional ^ 
and m have been scaled by a factor of M to match with the continuum limit. A little algebra shows that 7'* can be 
commuted through in the determinant above to yield the same overlap matrix of eq. ([7]) with exact chiral symmetry 
on the lattice. 
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IV. DISCUSSION 

The action in eq.([7]) leads to an overlap fermion determinant which is identical to that in the recent work (23| with 
fermionic sources in the overlap formalism of [24| . The main difference is, however, in the necessity of sources in (23j 
to define the chiral symmetry. Indeed, the chiral symmetry transformation there is local, defined as rotation of the 
sources while our eq. ([5]) is nonlocal, defined as the rotation of quark fields. The left-right symmetry is in-built in the 
formalism there whereas we needed to introduce the left-right projections in form of L- and i?-fields to do so. Our 
new fermion matrix is linear in /Lt^, similar to an earlier proposal by us p^ . This leads us to expect it to have the 
same good as well as bad properties. In particular, only its first derivative with respect to a/is is nonzero, all the rest 
being zero. As a consequence, the coefficients of the Taylor expansion of the baryonic susceptibility in /is, needed 
in estimating the location of the QCD critical point simplify considerably. On the other hand, the corresponding 
free theory itself has a /x^-divergence in the baryonic susceptibility in the continuum limit. It is easy to see that 



the successful prescription 25[ for local actions of introducing /zb as the fourth component of a constant Abelian 
gauge field to remove the divergence, will not work in this case since the nonlocal overlap Dirac matrix intertwines all 
four momentum components in a covariant manner. We have recently shown in case of the staggered fermions that 
the free theory divergence can be subtracted [1^ out successfully. In particular, the resultant ratios of the Taylor 
coefficients appear to be in good agreement with those where the subtraction is effected analytically by a change of 
action. While further investigations of the finite cut-off and finite volume effects are needed to establish it firmly, it 
may be hoped that a similar subtraction scheme will work for our above overlap quarks as well. Of course, it would be 
clearly desirable to modify even the action (O, without loss of its chiral symmetry, such that it has no /i^-divergences. 
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